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Variation of the fine structure constant in QSO spectra 
from coherent dark matter oscillations 
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ABSTRACT 

We consider the problem of the evolution of the fine structure coefficient a under the 
assumption that the scalar field coupling to the Maxwell term satisfies the condition 
rat 3> 1 for coherent dark matter oscillations. 

In this case we find that the coupling scale / in the leading order coupling 
—((/)/ A f)F flv F flv affects the cosmological evolution of a according to ln(a/ao) <x 
€(mpi/f) x ln(tanh(i/2r)/ tanh(t /2r)). A fit to the QSO observations by Murphy et 
al. yields / = £ X 2.12±g;|| x 10 5 m P; . Here m Pl = [SttGn)- 1 / 2 is the reduced Planck 
mass, and £ 2 = Q</>/Q m parametrizes the contribution of <j> to the matter density in the 
universe. 
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1 INTRODUCTION 

The question whether the value of Sommerfeld's fine struc- 
ture constant should actually be determined through the 
dynamics o f a sc alar field had been addressed already by 
iFierzl Jl95tjft and IJordanl (1959), whose investigations were 
partly motivated by Kaluza-Klein theory and by Dirac's 
proposal of a variability of constants over cosmological time 
scales. Nowadays it is well known that dynamical gauge 
couplings are predicted by string theory, and the rele- 
vant coupling e.g. of the heterotic string dilaton to gauge 
fields in four dimension s is of gravi tational strength / = 
(167rGjv) 1 /2 = y^m" 1 foicklll997ch . In an independent de- 
velopment Bckcnstcin ( 1982) had introduced a class of mod- 
els for dynamical a where the evolution of the fine structure 
constant is driven through couplings to energy densities. 

Therefore there was always theoretical interest in 
dynamical models for a, but in recent years Webb 
et al. also reported evidence for a variation of the 

time scales 



shift z and a = ceo = a(to). In iMurphv et all J2003I) re- 
sults are reported for 143 quasar apsorption lines, yielding 



x 10" 



1.75, but in our 



JWebb et alJ ll99Sl IMurphv et al.l 


20011: IWebb et alJ 1200 ll 


Murphy, Webb & Flambaum 


2003 


). The analysis of 128 


quasar absorption s 


ystems by 


Murphy. Webb & Flambaum 



( 2003) found an average increase 1 
— = ^±-l = -(0.543 ± 0.116) x HT 5 , 



(1) 



since redshift z = 1.67, i.e. over the last ~ 9.6 billion 
years. Here a z is the fine structure coefficient at red- 



We follow the standard sign convention Act = ct z — cxq, such 
that Aq < corresponds to an increase of a with time. 



Aa/a = -(0.57 ± 0.11) 
present analysis we used the well d ocumented sample from 
IMurphv. Webb fc Flambaum! fe003t . 

Dynamical gauge couplings can equivalently 
be expressed as dynamic al per meabilities, see e.g. 
Mag ueiio. Sandvik fc Kibble! d200ll) . Suppose q is a 
particular fixed value for the dynamical gauge coupling 
Q(x) (q will be further specified below). With the trans- 
formation qA /1 (x) = Q(x)An(x) the covariant derivatives 
can be written in terms of a variable or a constant gauge 
coupling 

D fl (x) = dp - iQ(x)Ap(x) = 9 M - iqA^x), 

but in the theory with the manifestly variable gauge cou- 
pling Q(x) the field strength tensor is 

Tv» = ^[D^,D V ] = ~0 M (QA,) - ^(QA^) 

while the field strength tensor with the constant coupling 
has the standard form F^ u — d^Av — d v A^ = (Q/q)^^. 
As a consequence the gauge theory with variable coupling 
and constant permeability appears as a gauge theory with 
constant coupling and variable permeability: 

-J-r^ix^ix) = — £ F^(x)F^(x). 

The dynamical coupling constant for charge Ze and the 
dynamical permeability in SI units are 



Q(x) = Ze(x)/h = 1Z\J ntoca(x)/h 
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and 



respectively. In the sequel we use units with h = c = 1. 

The scalar variable Q(x) may not have a canonically 
normalized kinetic term. Therefore a transformation Q{x) = 
Q((f>(x)) may be required if we want in leading order a stan- 
dard (d(j)) 2 term for the dynamics of Q: 



C 



4Q 2 W) 



^d^d^-V ■{<)>). 



(2) 



Here we assume that the potential V(4>) has a unique 
minimum V((j>o) at some value <j} = (j>o, and we parametrize 
the scalar field 4> such that <f>Q = 0. Furthermore, any non- 
vanishing term V^(O) would contribute to the cosmological 
constant and will not be considered as part of the energy 
density g^, stored in the scalar field 4>. The leading order 
expansion of the potential is then 



d 2 v{4>) 



(3) 



= 



With this proviso it seems prudent to choose q = Q(0) as the 
equilibrium value to which the gauge coupling should evolve 
due to the presence of the Hubble term 2>H4>. In leading order 
this implies the following parametrization for the coupling 
function Q((f>): 



^ 2 (4>) ~ q 2 ( I - j 



with the coupling scale denned accordingly 



1_ dQ 2 {4>) 

r,2 rIA 



(4) 



(•») 



0=0 



Besides the convention <j>o — for the equilibrium 
position this also implies a sign convention on the field 
4> if we require / > 0: is chosen as positive if it re- 
duces the fine structure cons tant in first order (see also e.g. 
iDamour fc Nordtvedd l|l993^ . and the observations of Mur- 
phy et al. then indicate that cj> is decaying from a positive 
value towards its equilibrium value (j>o — 0. 

Examples of specific coupling functions are provided e.g. 
by string theory or Kaluza-Klein theories: 



<2 2 0) = <? 2 exp(-0//), 



(6) 



and the Coulomb problem in these theories exh ibits an ultra - 
violet regularization at a scale 2 rj = q/(8nf) iDicklll997bl) . 

For the present investigation we will not specify the 
coupling function Q{4>) any further but only use the linear 
exp ansion @. 

lLandau fe Vucetichl i2002t) reconsidered the original 
Bekenstein model and concluded that it would not com- 
ply with the observ ations of Murphy et al. In a very in- 
teresting extension lOlive fc Pospelovl (I2002T) investigated 
a model where the dynamics of a massless scalar field 

2 Subsequently the abelian and non-abelian Coulomb prob- 
lem was also found to b e exactly solvable for other 



functions Q(d>) iDiclJ Il999t IChabab. Markazi fc SaidH |2000L 
ISlusarczvk fc Wereszczvnskil2001. 2003ft . and with mass couplings 
of d> iDicklll997dHBekensteirJl2002F 



<J> is mostly driven by its couplings to dark matter and 
the cosmological constant, and they analyzed compatibil- 
ity of Q with various constraints on variations of a. 
Our coupling paramete r / is related to the parameters 
M„ Cf and w in lOlive fc Pospelovl tOOtf) through / = 
M»/("f = V2luMpi/£f = and their results 

favor / > 10 3 A/p;, corresponding to a subgravitational cou- 
plin g strength of <j> to photons. 

iGardneJ J2003T) has recently discussed the implications 
of a mass term on the evolution of the fine-structure con- 
stant, and reported it to be consistent with mass val- 
ues for the scalar field <f> around m ~ Ho — I0~ 33 eV. 
Three crucial assumptions in Gardner's work are that the 
contribution of cf> to the dark matter density is negligi- 
ble, that / < mpi, and that < \(, m \m 2 Pl /f 2 < 10~ 5 , 
where £ m is the coupling of cfi to matter. A low mass 
value m ~ Hp was also prefer red in the recent work by 
lAnchordoaui fc Goldbergl i2005tl . who identified <j> with the 
quintessence field, and contrary to Gardner also assume d 
f > 10 3 Mpi in accordance with loiive fc Pospelovl feOOSft . 
ICopeland. Nunes fc Pospelovl teOO&t) also identified (f> with 
the quintessence and concluded that / ~ 10 5 mpi to fit the 
QSO data. However, this result did not comply with the 
Oklo constraint, and Copeland et al. proposed that a photon 
momentum dependence of / around 10 MeV may suppress 
the effects of dynamical a in nuclear reactions. 

In the present paper we propose yet another analysis of 
the implications of the results of Murphy et al. for a dy- 
namical fine structure constant. In particular we assume 
m > 10 -28 eV for the mass of the scalar field 4> generat- 
ing the evolution of a. Under this assumption a very weakly 
coupled field behaves like pressureless dust ever since dust 
domination, even though it may not satisfy the usual ther- 
mal dust condition 3 m > T. The virtue of m > 10" 28 eV 
for our present analysis is that under this condition we can 
use the late time behavior of 4>(t) for t 3> mT 1 to charac- 
terize the evolution of cf> ever since radiation-dust equality. 
Furthermore, our ignorance about evolution of (j> during ra- 
diation domination can be collected in a single parameter 4 
£ = [qqI ' Qm) 1 ^ 2 , and we perform a least squares fit of the 
ti me evolution of a in our mode l to th e a z values reported 
bv lMurphv. Webb fc Flambauml feOOSt . 

This explo res a completely different mass range than 
iGardnerl |2003). For mt 3> 1 the mass term generates tem- 
poral and spatial fluctuations of (j> at scales m _1 , but the 
Hubble expansion damps these oscillations oc t _1 , such that 
the amplitude of these oscillations is well below current ob- 
servational limits from laboratory based search experiments 
for variable a. Furthermore, with mt > 1 we will be able 
to use a virial theorem to eliminate m from the long term 
variation of <f>. The fit of the long term behavior of ait) de- 



3 This general result that non-thermal coherent oscillations be- 
have like cold d ark matter was observed for the first time 
in axion physics lAbbott fc Sikivielll983t iDine fc Fischlerlll983t 
iPreskill. Wise fc Wilczeklll983ft . 

4 ^ — 1 would imply that <fr is a dominant cold dark matter 
component. This possibilit y was pointed out for heavy dilatons 
(m S> T, dilaton wimps) bylGasperini fc Veneziand 11994) and by 

IDamour fc VilenkirJ il996ft. and for oscillations of light dilatons 
(T > m •» t' 1 ) bv lDickl <1997aft . 
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rived in Sec. [5] to the data of Murphy et al. then allows us 
to estimate the parameter //£. 

For cosmological parameters we rely on the recent eval- 
uation from the Wilkinson Microwave Anisotropy Probe 
(WMAP) and the Sloan Digital Sky Survey (SDSS) 
Ispereel et alJl2003l ; iTeemark etal]|2003^ . We use in partic- 
ular t he values from the "vanilla lite" model llTegmark et alJ 
2003): 



0.707™, ^ = 0.708™, 



(7) 



O a = 

to = 13.40lg:g Gyr. 

We included the errors given bv lTeemark et all J2f)03) for il- 
lustration, but do not use them for error propagation. They 
are negligible compared to the uncertainties in the a val- 
ues for the QSO absorption systems, which generate a la- 
uncertainty of about 22%, see Eq. 12011 below. 

The dependence of the WMAP CMB results on a is rel- 
atively weak in tha t it complies with 0.95 < ddec/cto < 1-02 
jRocha et alJl2003D . The dynamical evolution of a calcu- 
lated below implies that at the time of decoupling > 
(o-dec — a) /a > — 1 x 10 -4 , such that at this stage we can 
safely use WMAP results on cosmological parameters for the 
determination of //£. 

Scc.|2]recalls the relevant features of the dynamical evo- 
lution of a scalar field in an expanding universe and includes 
a virial theorem that will be useful in the analysis of dynam- 
ical gauge couplings. 

Our main result in Sec.[3]is an equation for the evolution 
of a(t) from the 1/1-7 coupling, and the fit to the results of 
Murphy et al. in Sec.0]yields //£. In Sec.|S]we will compare 
the time evolution of a in our model with the Oklo and 
meteorite constraints, and Sec. |S| contains our conclusions. 



2 THE COSMOLOGICAL EVOLUTION OF 4> 

With mt 3> 1 the mass parameter induces spatial and tem- 
poral fluctuations of <f>, and therefore of a. One might hope 
to use this to determine m from a Fourier decomposition of 
observations of a over cosmological distances. Our primary 
interest here is the coupling scale / of the scalar field to 
photons, and the strategy is to use a fit of the long term 
evolution of <j> in the expanding universe to the time varia- 
tion of a reported by Murphy et al. At this stage this allows 
us to infer a value for //£. 

<j> is usually assumed to have at most extremely weak 
matter couplings, and the long term evolution of very weakly 
coupled helicity states follows 



<f>(t) + 3H(t)4>(t) + m 2 </>(t) = 0, 



(8) 



with a corresponding evolution of the comoving energy den- 
3 (ft + m 2 ft)/2 



sity 
1 „ 



(3 ± 2 1 o ^ 
a <p + a m 



= ~^Ha a (m z 







(9) 



Note that in Eq. JHJ we used the convention 0o = 0, cf. ©. 
The solution of Eq. © for m — 0, 

dofti) oc a~ 3 (t) 

implies Q^it) oc a _6 (t), as appropriate for the ultrahard fluid 
component ge nerated by massless weakly coupled helicity 
states (see e.g. lDickl (120011) 1. 



The solution of Eq. © for m > and a ~ t 2/3 
<p(t) = (AJi (mt) + BYi(mt)j 
has asymptotics for mt 3> 1: 
4>(t) oc — cos(mt + (p), 

(>,/, oc t~ 2 oc a~ 3 . (10) 
This means that at late times ~ since 

just as for thermalized non-relativistic matter, but here even 
for m^T. 

Eq. 10 can be used to express the difference of the 
comoving kinetic and potential energy densities as a time 
derivative 

2a 3 H kln - 2a 3 H po t = a 3 ft - m 2 a 3 ft = ^(^H)- (11) 
This implies for the time limit 

H= lim - [ dtH(t) 
T Jo 

a virial theorem 



a 3 Hkin = a 3 TL P ot 



o 3 ^ = a z H = m 2 a 3 <t> 2 = a 3 ft. (12) 

However, note that at late times a 3 <fxj> oc t°, and therefore 
Eq. pip also yields 



ttkin — ttpoti 

i.e. 

|0| - m\<f>\, 

£></> ~ (p ■ 



(13) 

(14) 
(15) 



This relation can be used to trade the mass dependent \<j)\ for 
the energy density qj, in the late time evolution (t S> m _1 ) 
of a. 



3 THE COSMOLOGICAL EVOLUTION OF a 

From Eq. we have with a = ao = a(to): 

a(t) _ e 2 (ftt)) _ l 4> 
a e 2 f 

and with 1151 

(16) 

a / / 
We know from Eq. ffllUB and Eq. 1A5I in the Appendix that 

m ( a ° ^ 0*>o u\ Q<1>< A 
Q4,(t) = Q^, -7-T = Qm{t) = . 2 - . r - 

\a(t) J Q m ,o Qm.,0 sum (t/r) 

We denote the contribution from the dynamical gauge cou- 
pling to the matter density by £ 2 = f?</>,o/f?m,o, and find for 
the rate of change of the fine structure constant 



id(t)i _ i va 



a f sinh(t/r) 



(17) 
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Integration yields 

V«(*o) / ^ ^3/ \ v tanh(f /2r) / v ; 

The time constant is with A = 0.707f} c 

t = = 5.250 x 10 32 eV" 1 = 3.455 x 10 17 s = 10.95 Gyr. 

To match Eq. to the data from 

iMurphv. Webb fc Flambauml i2003h we also need the 
redshift-time relation for z <C 10 3 from 1A4H 



t z = T sinh 



i / sinh (t It) 



(19) 



lu - <,r> - <mU 1 fi+^i.5 °y r - 



(l + z) 1 - 5 
_! ^ 1.553 

k (l + *) 

Eqs. 1181 and 119H allow for a determination of the pa- 
rameter //£. 



4 THE COUPLING SCALE 

The parameter //£ was determined from a fit of 
the Eqs. 118^19^ to the a(z) values reported by 
IMurphv. W ebb & F lambauml J2003h . We set 



= In 



In 



a(tj) 
a (t ) 

/ tanh(t t /2r) 
I tanh(t /2r) 



and the minimal variance 

Ef Vi — sxj N " 
^ StH 

in the fit of 1181 1 occurs for a slope 
.2m P , J2i X iVi/ S yi 



3 C Vs/ E^?/^| 

with variance 

The errors in ln(a z /ao) are related to the errors in Aa/ao 
through 

SAa 

^ = — nr~ ■ 

ao + Aa 
This method yields a slope 

s = (5.443 ± 1.174) x 10~ 6 (20) 

corresponding to a coupling parameter 

j = 2.12±o;37 x 10 5 m Pi . (21) 

The variance per degree of freedom is Xdof ~ 1-067. 

The resulting asymptotic equilibrium value of the fine 
structure constant is 

lim a(t)/a(t ) = tanh(t /2r) _s ~ 1 + (3.3 ± 0.7) x 10~ 6 . 



The current rate of change of a from Eqs. 118I20H 
(4.0 ±0.9) x 10~ 17 yr _1 



a 4r sinh(to/T") 



(22) 



is within the bound s from current atomic clock experiments 
jMarion et al.ll2003D : 



-2.0 x 10~ 16 yr" 1 < - < 1.2 x 10~ 16 yr" 1 . 
a 



5 COMPARISON WITH THE OKLO AND 
METEORITE CONSTRAINTS 

With the cosmological parameters (J and Eq. 
11911 the closest QSO a bs orptio n system used in 



IMurphv. Webb fc Flambauml i2003fl corresponds to a 
distance of about 2.7 billion light years. A well known more 
recent constraint on variations of a over cosmological time 
scales comes from isotope abundances in the natural Oklo 
reactor, which had been active about 1.8 billion yea rs ago 
jShlvakhteJl97d iDamour fc pvson|ll996l : lFuiiill2003h . 
The recent evaluation bv lFuiiil (I2003i) yields a bound 



-(0.8 ± 1.0) x 10" 



(23) 



ctpkio — a 
a 

whereas insertion of 121H into I18H yields 

aokl ° - — = -(6.38±l-52) x 10" 7 . 
a 

The situation appears to be different with the 187 Re con- 
straints, which limit the evolution of a over the last 4 .6 Gyr 
jPeebles fc Dickelll962l : |PvsorJll97l lOlive et aLlliool : The 
most recent constraint from IOlivee^alTl2003T) is 

^^ = -(0.8±0.8)xl0- 6 , 
a 

while Eqs. 1211181 yield 



«4.6 



-(1.95 ±0.42) x 10" 



lOlive fc Pospelovl J20021) and iGardnerl (I2003T) could 
fit their externally driven models for the evolution 
of a to both the QSO data and the O kl o con - 
straint, whereas [A nchord oaui fc Goldberd i2003ft . 
ICopeland. Nunes fc Pos pelovl j2003t) a nd w e find a lower 
value of a. However. iMotafcBarrowl l)2003|) have recently 
pointed out that the local variation of a in virialized over- 
densities like our own can relax the Oklo constraint by a 
factor 10-100, because on the one hand a dynamical a would 
be expected to have a higher value in overdensities, while on 
the other hand virialization of overdensities slows down the 
evolution of a. The qualitative picture emerging from this 
is that in overdensities a evolves from a higher initial value 
after virialization, but at slower pace, whence it approaches 
again the value in the low-density background universe. 
This can explain discrepancies between astrophysical and 
geochemical observations. 



6 CONCLUSIONS 

We have found the equation l|18p for the dynamical time evo- 
lution of a due to the coupling -~(4>/4 : f)F fJ ' v _F M „ of electro- 
magnetic fields to a very weakly coupled massive scalar field 
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with mt 3> 1. A fit of this equation to the quasar absorp- 
tion data reported by iMurphv. Webb fc Flambauml J2003f) 
yieids the value 12111 . where £ = yQ^o/Omfi parametrizes 
the contribution of the scalar field <f) to the matter density. 

Within this model the evolution of a reported by Mur- 
phy et al. appears to be slow due to a small coefficient 
s = (5.443±1.174)xlCn 6 in Eq. USJl: The fine-structure con- 
stant varied so little since z — 3.66 because the <f> abundance 
is small and the ^-7 coupling is very weak and presumably 
of sub gravitational stren g th, in agreement with the anal- 
yses o f lOliye fc Posoelovl <2002l): lAnchordoaui &; Goldberg 
(l2003h and lCooeland. Nunes fc Pospelovl <2003h . 

Our coherent oscillation model for dy- 
namical a (and other self-driven m odels of 

dynamical a iAnchordooui fc Goldberg 1 l2003t 

ICopeland. Nunes fc Pospelovl 12003^ 7 still seems to predict 
a too small value of a at the time when the Oklo n atural 
reactor was active. However, as lMota fc Barrowl J2003T) have 
pointed out, at low redshift predictions of faster evolution of 
q from astrophysical observations are to be expected due to 
spatial variations in the presence of local overdensities. Note 
that this dues not affect the nu merical results 120J21|22]I 
(or th e cor responding results of jA^ichordocjui fc Goldberg 
120031) and ICopeland. Nunes fc PospeloXil2o'o3h '). since the 
resulting discrepancy of qo on Earth and far away from 
virialized objects is smaller than the number of significant 
figures reported 5 . 
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APPENDIX A: THE SCALE FACTOR IN THE 
ACDM UNIVERSE 

The evolution of the scale factor in a spatially flat ACDM 
universe follows from direct integration of the corresponding 
Friedmann equation 



h_ _ Qm + A 

a 2 3m 2 pl 
after insertion of 

Qm(t) = Qm.Q 

Integration from to to t yields 



(Al) 



o(t) 



-(t-to) = In 



2mPl \ V^<4 + V( A + ern,o)d% / 

Solving for the scale factor yields 



H sinh^LJ"^ (A2) 



.fflo 

Qm,0 

A + ^/A 2 + Kg m ,o 
with the time constant 

r=*™ (A3) 
V3A 

We can simplify our result 1A2I because the highest 
redshift z = 3.66 used in the analysis still corresponds to 
an age t(z) ~ 1.7 Gyr 2> t Bq much larger than the time 
t eq ~ 1.3 x 10 5 yr of matter radiation equality. At times 
3> t e q the modification of the time evolution of the scale 
factor during the very early radiation dominated era can be 
neglected, and one can integrate Eq. iAll from ti = and 
still get an extremely good approximation for t S> t eq . This 
yields 

/ \ 2/3 

a _ ( sinh(t/r) \ . 



ao y sinh(io/i") 
and 

Qm = ^JWrj- (A5) 
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